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We study nuclear effects of charged current deep inelastic neutrino-iron scattering in the framework of a 
analysis of parton distribution functions (PDE?). We extract a set of iron PDI% which are used to compute 
a;sj-dependent and Q^-dependent nuclear correction factors for iron structure functions which are required in 
global analyses of free nucleon PDFs. We compare our results with nuclear correction factors from neutrino- 
nucleus scattering models and correction factors for ^^-iron scattering. We find that, except for very high XBj, 
our correction factors differ in both shape and magnitude from the correction factors of the models and charged- 
lepton scattering. 



1. Impact of Nuclear Corrections on PDFs 

The high statistics measurements of neutrino 
deeply inelastic scattering (DIS) on heavy nuclear 
targets has generated significant interest in the 
literature since these measurements provide valu- 
able information for global fits of parton distribu- 
tion functions (PDF^). It is necessary to use both 
Charged Current (CC) probes and Neutral 
Current (NC) {7, Z} probes to disentangle the 
separate PDF flavor components. Toward this 
goal, the use of nuclear targets is unavoidable due 
to the weak nature of the {W*, Z} interactions, 
and this complicates the extraction of free nu- 
cleon PDFs because model-dependent corrections 
must be applied to the data. 

In early PDF analyses, the nuclear corrections 
were static correction factors without any (sig- 
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nificant) dependence on the energy scale Q, the 
atomic number A, or the specific observable. The 
increasing precision of both the experimental data 
and the extracted PDFs demand that the applied 
nuclear correction factors be equally precise as 
these contributions play a crucial role in deter- 
mining the PDFs. 

In this study we reexamine the source and size 
of the nuclear corrections that enter the PDF 
global analysis, and quantify the associated un- 
certainty. Additionally, we provide the founda- 
tion for including the nuclear correction factors 
as a dynamic component of the global analysis so 
that the full correlations between the heavy and 
light target data can be exploited. 

A recent study[T] analyzed the impact of new 
data sets from the NuTeV, Chorus, and E-866 
Collaborations on the PDFs. This study found 
that the NuTeV data set (together with the model 
used for the nuclear corrections) pulled against 
several of the other data sets, notably the E-866, 
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BCDMS and NMC sets. Reducing the nuclear 
corrections at large values of x reduced the sever- 
ity of this pull and resulted in improved values. 
These results suggest on a purely phenomenolog- 
ical level that the appropriate nuclear corrections 
for I'-DIS may well be smaller than assumed. 

2. Global Analysis Framework 

To investigate this question further, we use the 
high-statistics zz-DIS experiments to perform a 
dedicated PDF fit to neutrino-iron data.|2] Since 
we first will study iron alone and will not (ini- 
tially) combine the data with measurements on 
different target materials, we need not make any 
assumptions about the nuclear corrections; this 
side-steps a number of difficulties. |3] [1] 3] While 
this approach has the advantage that we do not 
need to model the >l-dependence, it has the draw- 
back that the data from just one experiment will 
not be sufficient to constrain all the parton dis- 
tributions; therefore, other assumptions must en- 
ter the analysis. The theoretical framework will 
roughly follow the CTEQ6 analysis of free pro- 
ton PDFs.[5] We outline the key features of our 
analysis below, and focus on the issues specific to 
our study of NuTeV neutrino-iron data in terms 
of nuclear parton distribution functions. 

2.1. Basic formalism 

For our PDF analysis, we will use the general 
features of the QCD-improved parton model and 
the analyses as outlined in Ref. [5]. We adopt 
the framework of the recent CTEQ6 analysis of 
proton PDFs where the input distributions at the 
scale Qo = 1-3 GeV are parameterized as 



As 



xh{x, Qo) = Aox^ni - a;)-^^e'^^"(l + e^'^x) 

for i — {uv, dv,g, u + d,s, s}, and 

xMx, Qo) - Aox^' (1 - x)^' + (1 + A3x)il - x)^^ 

for i = {d/u} where Uy and are the up- and 
down-quark valence distributions, u, d, s, s are 
the up, down, strange and anti-strange sea dis- 
tributions, and g is the gluon. Furthermore, 
the fi — f^^^ denote parton distributions of 



hound protons in the nucleus A, and the vari- 
able < a; < A is defined as x := Axa where 
xa = /"^Pa • 9 is the usual Bjorken variable 
formed out of the four-momenta of the nucleus 
and the exchanged boson. This parameterization 
is designed for < x < 1 and we here negleclj^the 
distributions at a; > 1. Note that the condition 
/i(a; > 1, Q) = is preserved by the DGLAP evo- 
lution and has the effect that the evolution equa- 
tions and sum rules for the /f are the same as 
in the free proton case. 

The PDFs for a nucleus (A, Z) are constructed 

as 



/^(x,Q)=^/f/^(x,Q)-H 



where we relate the distributions inside a bound 
neutron, f^^^{x,Q), to the ones in a proton by 
assuming isospin symmetry. The nuclear struc- 
ture functions are given by parallel relations such 
that they can be computed in next-to-leading or- 
der as convolutions of the nuclear PDFs with the 
conventional Wilson coefficients, i.e., generically 

In order to take into account heavy quark 
mass effects we calculate the relevant structure 
functions in the ACOT schemedH [12] in NLO 
QCD.[T2 Finally, the differential cross section for 
charged current (anti-)neutrino-nucleus scatter- 
ing is given in terms of three structure functions: 



d^a 
dx dy 
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where the ('— ') sign refers to neutrino (anti- 
neutrino) scattering and where Gp is the Fermi 
constant, M the nucleon mass, and E the energy 
of the incoming lepton (in the laboratory frame). 

2.2. Methodology 

The basic formalism described in the previ- 
ous sections is implemented in a global PDF fit- 

^ While the nuclear PDFs can be finite for x > 1, the 
magnitude of the PDFs in this region is negligible for the 
purposes of the present study (c/., Refs. [6ll71[8l[9l[T0]1. 
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ting package, but with the difference that no nu- 
clear corrections are applied to the analyzed data; 
hence, the resulting PDFs are for a bound proton 
in an iron nucleus. The parameterization pro- 
vides enough flexibility to describe current data 
sets entering a global analysis of free nucleon 
PDFs; given that the nuclear modifications of the 
x-shape appearing in this analysis are modest, 
this parameterization will also accommodate the 
iron PDFs. 

Because the neutrino data alone do not have 
the power to constrain all of the PDF compo- 
nents, we will need to impose some minimal set 
of external constraints. For example, our results 
are rather insensitive to the details of the gluon 
distribution with respect to both the overall 
and also the effect on the quark distributions. 
The nuclear gluon distribution is very weakly con- 
strained by present data, and a gluon PDF with 
small nuclear modifications has been found in the 
NLO analysis of Ref. We have therefore fixed 
the gluon input parameters to their free nucleon 
values. For the same reasons the gluon is not sen- 
sitive to this analysis, fixing the gluon will have 
minimal effect on our results. Furthermore, we 
have set the d/u ratio to the free nucleon result 
assuming that the nuclear modifications to the 
down and up sea are similar such that they can- 
cel in the ratio. This assumption is supported by 
Fig. 6 in Ref. [Ti] . 



3. Analysis of iron data 

3.1. Iron Data Sets 

We determine iron PDFs using the recent 
NuTeV differential neutrino and anti-neutrino 
DIS cross section data.jlS] In addition, we include 
NuTeV/CCFR dimuon data |16] which are sensi- 
tive to the strange quark content of the nucleon. 
There are other measurements of neutrino-iron 
DIS available in the literature from the CCFR 
Ha Ha Sg |2Q], CDHS[21] and CDHSW|22] col- 
laborations; see, e.g., Ref. [23] for a review. There 
is also a wealth of charged lepton-iron DIS data 
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Figure 1. Parton distributions for iron. The cen- 
tral PDF from fit 'A2' is shown by the solid line. 
The dashed lines depict parton distributions con- 
structed with ^ = 56 and Z = 26 using the 
Base-1 free-proton PDI^. Additional results are 
shown from HKN04 0, (NLO) HKN07 |8], and 
(DS) [11]. The vertical line marks the lower limit 
of the data in the x variable. 



including SLAC[23 and EMC [25] [26] |^ For the 
initial study we limit our analysis to the NuTeV 
experiment alone; we will compare and contrast 



*C/. the Durham HEP Databases for a complete listing: 
http : //www- spires . dur . ac .uk/hepdata/ 
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different experiments in a subsequent study. 

3.1.1. PDF Reference Sets 

For the purposes of this study, we use two dif- 
ferent reference sets of free-proton PDFs which 
we denote 'Base-1' and 'Base-2'. 

Since we focus on iron PDFs and the associated 
nuclear corrections, we need a base set of PDFs 
which are essentially free of any nuclear effects; 
this is the purpose of the Base-1 reference set [1]. 
Therefore, to extract the Base-1 PDFs we omit 
the CCFR and NuTeV data from our fit so that 
our base PDFs do not contain any large residual 
nuclear corrections]!] The absence of such nuclear 
effects will be important when we extract the nu- 
clear corrections factors. 

The Base-2 PDFs are essentially the 
CTEQ6.1M PDFs with a modified strange PDF 
introduced to accommodate the NuTeV dimuon 
data]^ In the manner of the CTEQ6.1M PDF's, 
the Base-2 fit does not apply any deuteron correc- 
tions to the data; this is in contrast to the Base-1 
PDFs. Also, the Base-2 fit does include the CCFR 
data that has been corrected to a free nucleon 
using charged-lepton correction factors. [20] 

By comparing the free-proton PDF 'Base-1' 
and 'Base-2' sets with the iron PDF sets, we can 
gauge the size of the nuclear effects. Furthermore, 
differences between observables using the 'Base-1' 
respectively the 'Base-2' reference sets will indi- 
cate the uncertainty due to the choice of the free- 
proton PDF. 

3.1.2. Comparison of the Fits with Data 

Specifically, we determine iron PDFs using the 
recent NuTeV differential neutrino (1371 data 

^We do retain the deuteron data as this has only a small 
correction over the central a;-range,[3l [l]. The deuteron 
correction has been applied in the Base-1 fit. Also, for 
the Drell-Yan Cu data (E605), the expected nuclear cor- 
rections in this kinematic range are small (a few percent) 
compared to the overall normalization uncertainty (15%) 
and systematic error (10%). 

®These PDFs have been determined from a fit to the same 
data set as in the CTEQ6 analysis with the addition of 
the NuTeV dimuon data. The changes to the strange sea 
induce only minor changes to the other fit parameters; this 
has a minimal effect on the particular observables (da, F2) 
we examine in the present study. 



points) and anti-neutrino (1146 data points) 
DIS cross section data,jT5] and we include 
NuTeV/CCFR dimuon data (174 points) which 
are sensitive to the strange quark content of the 
nucleon. Using the ACOT scheme, we impose 
kinematic cuts of Q > 2 GeV and W > 3.5 GeV, 
and obtain a good fit with a oi 1.35 per data 
point; we identify this fit as 'A2.'[2] 

3.2. Iron PDFs 

We now examine the nuclear (iron) parton dis- 
tributions f.f{x,Q^) in Figure m which shows the 
PDFs from fit 'A2' at our input scale Qq = rric = 
1.3 GeV versus x. For an almost isoscalar nu- 
cleus like iron the u and d distributions are very 
similar. Therefore, we only show the it„ and u 
partons, together with the strange sea]^ As ex- 
plained above, the gluon distribution is very sim- 
ilar to the familiar CTEQ6M gluon at the input 
scale such that we don't show it here. In order 
to indicate the constraining power of the NuTeV 
data, the band of reasonable fits is depicted. The 
fits in this band were obtained (as outlined above) 
by varying the initial conditions and the num- 
ber of free parameters to fully explore the solu- 
tion space. All the fits shown in the band have 

/DOF within 0.02, which roughly corresponds 
to a range of Ax^ ~ 50 for the 2691 data points. 

As can be seen in Figure [l] the Uy distribu- 
tion has a very narrow band across the entire x- 
range. The up- and strange-sea distributions are 
less precisely determined. At values of x down to, 
say, X ~ 0.07 the bands are still reasonably well 
confined; however, they open up widely in the 
small- cc region. Cases where the strange quark 
sea lies above the up-quark sea are unrealistic, 
but are present in some of the fits since this region 
{x < 0.02) is not constrained by data. We have 
included the curves for our free-proton Base-1 
PDFs (dashed), as well as the HKN04 |7] (dot- 
ted), the NLO HKN07 g] (dotted-dashed), and 

^While iron is roughly isoscalar, other nuclear PDFs 
can exhibit larger differences between the u and d 
distributions — the extreme case being the free-proton 
PDF. When comparing PDFs we must keep in mind that 
it is ultimately the structure functions which are the phys- 
ical observables. 
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DS [H] (dot-dashed) nuclear PDFsj^ 

The comparison with the Base-1 PDFs is 
straightforward since the same theoretical frame- 
work (input scale, functional form, NLO evolu- 
tion) has been utilized for their determination. 
Therefore, the differences between the solid band 
('A2') and the dashed line (Base-1) exhibit the 
nuclear effects, keeping in mind that the free- 
proton PDFs themselves have uncertainties. 

For the comparison with the HKN04 distribu- 
tions, it should be noted that a SU(3)-fiavor sym- 
metric sea has been used; therefore, the HKN04 
strange quark distribution is larger, and the light 
quark sea smaller, than their Base-1 PDF coun- 
terparts over a wide range in x. Furthermore, the 
HKN04 PDFs are evolved at leading order. 

In a recent analysis, the HKN group has pub- 
lished a new set of NPDFs (HKN07) including un- 
certainties |8]. They provide both LO and NLO 
sets of PDFs, and we display the NLO set. These 
PDFs also use a more general set of sea distribu- 
tions such that u{x) d{x) ^ s(a;) in general. 

The DS PDFs are Hnked to the GRV98 PDFs 
|27] with a rather small radiatively generated 
strange sea distribution. Consequently, the light 
quark sea is enhanced compared to the other sets. 
Additionally, the DS sets are evolved in a 3-fixed- 
fiavor scheme in which no charm parton is in- 
cluded in the evolution. However, at the scale 
Q = rric of Fig. [ijthis is of no importance. 

4. Nuclear Correction Factors 

In the previous section we analyzed charged 
current ^-Fe data with the goal of extracting the 
iron nuclear parton distribution functions. In this 
section, we now compare our iron PDFs with the 
free-proton PDFs (appropriately scaled) to infer 
the proper heavy target correction which should 
be applied to relate these quantities. 

Within the parton model, a nuclear correction 
factor R[0] for an observable O can be defined as 



follows: 



*In a recent publication, Eskola et al. [To] perform a global 
reanalysis of their ESK98 [9] nuclear PDFs. While we do 
not present a comparison here, the results are compatible 
with those distributions displayed in Fig^H a comparison 
can be found in Figs. 10 and 11 of Ref. [tO". 



R[0] = 



©[NPDF] 



(1) 



0[free] 

where ^^[NPDF] represents the observable com- 
puted with nuclear PDFs, and ^^[free] is the same 
observable constructed out of the free nucleon 
PDFs. Clearly, R can depend on the observable 
under consideration simply because different ob- 
servables may be sensitive to different combina- 
tions of PDFs. 

This means that the nuclear correction factor 
R for and F^^ will, in general, be different. 
Additionally, the nuclear correction factor for F2 
will yield different results for the charged current 
v-Fe process {W^ exchange) as compared with 
the neutral current £^-Fe process (7 exchange). 
Schematically, we can write the nuclear correction 
for the DIS structure function F2 in a charged 
current (CC) v-A process asj^ 



R'cciF2;x,Q^) 



(2) 



and contrast this with the neutral current (NC) 
l^-A process: 

i?^^(F2;a:,g2)~ 
(-i)^[rf^-fJ^ + ...]+(+i)^K + ^^ + ...] 

(-i)' [d0 + rf0 + ...] + (+§)' [«0 + 00 + ...] ' 

where the superscript " 0" denotes the "free nu- 
cleon" PDF which is constructed via the relation: 

//(^,Q)-|/f(x,Q) + ^^^^/r(2^,Q) -(3) 

Clearly, the i?-factors depend on both the kine- 
matic variables and the factorization scale. Fi- 
nally, we note that Eq. ([ij is subject to uncertain- 
ties of both the numerator and the denominator. 

We will now evaluate the nuclear correction fac- 
tors for our extracted PDF^, and compare these 
with selected results from the literature l28ll29lFn 



^The corresponding anti-neutrino process is obtained with 
a « <-» d interchange. 

^"Note that our comparison with the Kulagin-Petti model 
is based on the work in Ref l28l. 
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Figure 2. Parameterization for the neutral cur- 
rent charged lepton structure function F^'^/F^. 
For comparison we show experimental results 
from the BCDMS collaboration (BCDMS-85 [M], 
BCDMS-87 |3T]) and from experiments at SLAC 
(SLAC-E049 |3l], SLAC-E139 0, and SLAC- 
E140 [21] )• Normalization uncertainties of the 
data have not been included. 
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Figure 3. Nuclear correction factor R for the dif- 
ferential cross section cPa/dx dQ^ in charged cur- 
rent neutrino- Fe scattering at = 5 GeV^. Re- 
sults are shown for the charged current neutrino 
(solid lines) and anti-neutrino (dashed lines) scat- 
tering from iron. The upper (lower) pair of curves 
shows the result of our analysis with the Base-2 
(Base-1) free-proton PDFs. 



Because we have extracted the iron PDFs from 
only iron data, we do not assume any particu- 
lar form for the nuclear yl-dependence; hence the 
extracted R[0] ratio is essentially model indepen- 
dent. 

4.1. F^^/Ff NC charged lepton scattering 

We will also find it instructive to compare our 
results with the as extracted in neu- 

tral current charged-lepton scattering, i^-Fe. In 
Fig. [2] we compare the experimental results for 
the structure function ratio for the fol- 

lowing experiments: BCDMS-85 H, BCDMS-87 
EH, SLAC-E049 [32], SLAC-E139 |3], SLAC-140 
|2i] . The curve (labeled SLAC/NMC parame- 
terization) is a fit to this data. While there is 
a spread in the individual data points, the pa- 
rameterization describes the bulk of the data at 
the level of a few percent or better. It is im- 
portant to note that this parameterization is in- 
dependent of atomic number A and the energy 
scale |33]; this is in contrast to the results we 
will derive using the PDFs extracted from the nu- 



clear dataj^ Additionally, we note that while this 
parameterization has been extracted using ratios 
of F2 structure functions, it is often applied to 
other observables such as Fi,3,l or da. We can 
use this parameterization as a guide to judge the 
approximate correspondence between this neutral 
current (NC) charged lepton DIS data and our 
charged current (CC) neutrino DIS data. 

4.2. R Factors for d'^a/dxdQ'^ 

We begin by computing the nuclear correction 
factor R according to Eq. ([ij for the neutrino dif- 
ferential cross section as this represents the bulk 
of the NuTeV data included in our fit. More pre- 
cisely, we show i?-factors for the charged current 
cross sections d^a/dxdQ^ at fixed Q^. Our re- 
sults are displayed in Fig. |3]for = 5 GeV^ and 
a neutrino energy Ei, — 150 GeV which implies, 
due to the relation — 2MEi,xy, a minimal x- 
value of Xmin — 0.018. The solid (dashed) lines 

^^In particular, we will find for large x (> 0.5) and Q 
comparable to the proton mass the target mass corrections 
for F^'^/F^ are essential for reproducing the features of 
the data; hence the Q dependence plays a fundamental 
role. 
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correspond to neutrino (anti-neutrino) scattering 
using the iron PDFs from the 'A2' fit. 

We have computed R using both the Base-1 
and Base-2 PDFs for the denominator of Eq. ([ij; 
recall that Base-1 includes a deuteron correction 
while Base-2 uses the CCFR data and does not 
include a deuteron correction. The difference be- 
tween the Base-1 and Base-2 curves is approxi- 
mately 2% at small x and grows to 5% at larger 
X, with Base-2 above the Base-1 results. The dif- 
ference of these curves, in part, reflects the uncer- 
tainty introduced by the proton PDF. [2] As this 
behavior is typical, in the following plots (Figs.|4j 
we will only show the Base-1 results. We also ob- 
serve that the neutrino (anti-neutrino) results co- 
incide in the region of large x where the valence 
PDFs are dominant, but differ by a few percent 
at small x due to the differing strange and charm 
distributions. 

4.3. R Factors for F^{x,Q^) and F^{x,Q^) 

We now compute the nuclear correction fac- 
tors for charged current neutrino-iron scatter- 
ing. The results for I'-Fe and those of D-Fe are 
shown in Fig. |4] The numerator in Eq. ([iJ has 
been computed using the nuclear PDF from fit 
'A2', and for the denominator we have used the 
Base-1 PDFs. For comparison we also show the 
correction factor from the Kulagin-Petti model 
[281 [29] (dashed-dotted), and the SLAC/NMC 
curve (dashed) which has been obtained from an 
A and Q^-independent parameterization of cal- 
cium and iron charged-lepton DIS data. 

Due to the neutron excess in iron, both our 
curves and the KP curves differ when comparing 
scattering for neutrinos and anti-neutrinos; the 
SLAC /NMC parameterization is the same in both 
figures. For our results (solid lines), the difference 
between the neutrino and anti-neutrino results is 
relatively small, of order 3% at a; = 0.6. Con- 
versely, for the KP model (dashed-dotted lines) 
the I'-D difference reaches 10% at a; ~ 0.7, and 
remains sizable at lower values of x. 

To demonstrate the dependence of the R fac- 
tor on the kinematic variables, in Fig. |4]we have 
plotted the nuclear correction factor for two sep- 



arate values of Q^. Again, our curves and the 
KP model yield different results for different 
values, in contrast to the SLAC/NMC parame- 
terization. 

Comparing the nuclear correction factors for 
the F2 structure function with those obtained for 
the differential cross section (Fig. |3}, we see these 
are quite different, particularly at small x. Again, 
this is because the cross section (Pa is comprised 
of a different combination of PDFs than the F2 
structure function. In general, our i?-values for 
F2 lie below those of the corresponding i?- values 
for the cross section da at small x. Since da is 
a linear combination of F2 and F3, the i?- values 
for ^3 (not shown) therefore lie above those of 
F2 and da. Again, we emphasize that it is im- 
portant to use an appropriate nuclear correction 
factor which is matched to the particular observ- 
able. 

As we observed in the previous section, our re- 
sults have general features in common with the 
KP model and the SLAC/NMC parameteriza- 
tion, but the magnitude of the effects and the 
a:;- region where they apply are quite different. 
Our results are noticeably flatter than the KP 
and SLAC/NMC curves, especially at moderate- 
X where the differences are significant. The gen- 
eral trend we see when examining these nuclear 
correction factors is that the anti-shadowing re- 
gion is shifted to smaller x values and any turn- 
over at low X is minimal given the PDF uncer- 
tainties. In general, these plots suggest that the 
size of the nuclear corrections extracted from the 
NuTeV data are smaller than those obtained from 
charged lepton scattering (SLAC/NMC) or from 
the set of data used in the KP model. We will 
investigate this difference further in the following 
section. 

4.4. F^^/F^ from Iron PDFs 

Since the SLAC/NMC parameterization was fit 
to F2^ /F2 for charged-lepton DIS data, we can 
perform a more balanced comparison by using our 
iron PDFs to compute this same quantity. The 
results are shown in Fig. [5] where we have used 
our iron PDFs to compute -Flf^'^, and the Base-1 
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Figure 4. Nuclear correction factor R for the structure function F2 in neutrino and anti-neutrino scat- 
tering from Fe for = {5,20}GeV^. The soUd curve shows the result of our analysis of NuTeV data; 
the uncertainty from the fit is represented by the shaded (yellow) band. For comparison we show the 
correction factor from the Kulagin-Petti model (dashed-dot line), [28] HKN07 (dashed-dotted line),|8] 
and the SLAC/NMC parametrization (dashed line). 



and Base-2 PDFs to compute . 

As with the nuclear correction factor results 
of the previous section, we find our results have 
some gross features in common while on a more 
refined level the magnitude of the nuclear cor- 
rections extracted from the CC iron data differs 
from the charged lepton data. In particular, we 
note that the so-called "anti-shadowing" enhance- 
ment aX X ^ [0.06 — 0.3] is not reproduced by the 
charged current (anti-) neutrino data. Examin- 
ing our results among all the various R\0] cal- 
culations, we generally find that any nuclear en- 
hancement in the small x region is reduced and 
shifted to a lower x range as compared with the 
SLAC/NMC parameterization. Furthermore, in 



the limit of large x {x >^ 0.6) our results are 
slightly higher than the data, including the very 
precise SLAC-E139 points; however, the large the- 
oretical uncertainties on in this x-region make 
it difficult to extract firm conclusions. 

This discussion raises the more general ques- 
tion as to whether the charged current {v-Fe) 
and neutral current (£*-Fe) correction factors 
are entirely compatible [TSl l3i] l35] l36l l371 l38]. 
There is a priori no requirement that these be 
equal; in fact, given that the v-Fe process in- 
volves the exchange of a and the i^-Fe pro- 
cess involves the exchange of a 7 we necessarily 
expect this will lead to differences at some level p| 

^■^In Ref. f39l, Brodsky and collaborators posit a non- 



Nuclear PDFs 
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Figure 5. Predictions (solid and dashed line) for 
the structure function ratio FJ^'^/F^ using the 
iron PDFs extracted from fits to NuTeV neutrino 
and anti-neutrino data. The SLAC/NMC pa- 
rameterization is shown with the dot-dashed line. 
The structure function F^ in the denominator 
has been computed using either the Base-2 (solid 
line) or the Base-1 (dashed line) PDFfe. 



4.5. Future Studies 

It is important to resolve whether the differ- 
ences we observe in Fig. [5] arise from the un- 
certainty of the nuclear corrections, or if they 
are genuinely a consequence of NC/CC effects. 
A combined analysis of CC neutrino and NC 
charged-lepton data sets will shed more light on 
these issues. To best address these questions, we 
need to include the nuclear dimension (parame- 
terized by the nuclear A value) as a dynamic com- 
ponent of the global fit; this will allow us to fit 
both the CC V7*-exchange processes at large A, 
as well as the NC 7-exchange processes at small 
A in a coherent framework. Figure |6] presents 
an illustrative example of how the PDFs can be 
extended to incorporate the necessary A depen- 
dence to implement such a program. This ex- 
tended analysis with additional data sets is in 
progress, and should help clarify these questions. 



universal nuclear anti-shadowing mechanism which yields 
different effects for CC and NC scattering. 




0.02 0.05 0.10 0.20 0.50 1.00 



Figure 6. Illustration of the gluon PDF x ff{x, fi) 
vs. a; as a function of the nuclear A. 



5. CONCLUSIONS 

While the nuclear corrections extracted from 
charged current f-Fe scattering have similar 
characteristics as the neutral current l^-Fe 
charged-lepton results, the detailed x and be- 
havior is quite different. This observation raises 
the deeper question as to whether the charged 
current and neutral current nuclear correction 
factors may be substantially different. This 
present study of the iron PDFs provides a founda- 
tion for a general investigation (involving a vari- 
able A parameter) that can definitively address 
this topic. Resolving these questions is essential 
if we are to reliably use the plethora of nuclear 
data to obtaining free-proton PDFs which form 
the basis of the LHC analyses. 
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